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We present results of a systematic investigation of the evolution of low-mass (0.35 M Q , 0.40 M Q and 
0.65 M Q ) helium donors in semidetached binaries with accretors - white dwarfs. In the initial models 
of evolutionary sequences abundance of helium in the center 0. 1 ~ Y c < 0.98. Results of computations 
may be applied to the study of the origin and evolutionary state of AM CVn stars. It is shown that 
the minimum orbital periods of the systems only weakly depend on the total mass of the system and 
evolutionary state of the donor at RLOF and are equal to 9-1 1 min. The scatter in the mass-exchange 
rates at given P orb in the range P m b,min < Porb ~ 40 min. does not exceed ~ 2.5. At P 0lb > 20 min 
mass-losing stars are weakly degenerate homogeneous cooling objects and abundances of He, C, N, O, 
Ne in the matter lost by them depends on the extent of He-depletion at RLOF. For the systems which 
are currently considered as the most probable model candidates for AM CVn stars with helium donors 
these abundances are Y > 0.4, X c ~ 0.3, X < 0.25, X N < 0.5 x 10" 2 . At P orb > 40 min. the timescale 
of mass-loss begins to exceed thermal time-scale of the donors, the latter begin to contract, they 
become more degenerate and, apparently, "white-dwarf" and "helium-star" populations of AM CVn 
stars merge. 
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1 Introduction 



Non-degenerate helium stars in close binaries (CB) form in so-called "case B" of mass-exchange, 
when stars with mass M > (2.3 - 2.5) M© overflow their Roche lobes in the hydrogen- shell burning 
stage (Kippenhahn and Weigert, 1967; Paczyhski, 1967b). Masses of helium stars are > 0.32 M Q 
(Iben, 1990; Han et al., 2002). 

For the current notions on stellar evolution, importance of low-mass helium stars is defined by 
the possibility of formation of pairs containing white dwarfs or neutron stars accompanied by helium 
stars in the course of evolution of CB. Orbital periods of such systems may be so short, that angular 
momentum loss (AML) via gravitational wave radiation (GWR) enables Roche-lobe overflow (RLOF) 
by helium star before helium exhaustion in the core of the latter and, under proper conditions, stable 
mass-transfer is possible (Savonije et al., 1986; Iben and Tutukov, 1987; see also Yungelson, 2005a 
and Postnov and Yungelson, 2006). 

Semidetached pairs of neutron stars with helium-star companions were not observed as yet, 
though, the first computation of the evolution of a CB with a non-degenerate helium donor was per- 
formed just for such a pair (Savonije et al. 1986). On the other hand, it is assumed that in the case 
when in a semidetached binary helium star is accompanied by a white dwarf, such a system may be 
observed as an ultra-short period cataclysmic variable of the AM CVn type (Savonije et al., 1986). 
Other hypothetical scenarios for formation of AM CVn stars assume stable mass-exchange between 
white dwarfs (Paczyhski, 1967a) or mass-loss by a remnant of low-mass (< 1.5 M Q ) main-sequence 
star which overflowed its Roche lobe after exhaustion of a significant fraction of the hydrogen in its 
core (X c < 0.4) or even immediately after formation of a low-mass (~ 0.01 M Q ) helium core (Tutukov 
et al., 1985; Podsiadlowski et al., 20030 

AM CVn stars are important for the physics in general, since it is expected that, thanks to their 
extremely short orbital periods, they, along to the close detached short-period white dwarf pairs, will 
be among the first objects which will be able to detect space-born gravitational wave antenna LISA 
(Evans et. al., 1987; Nelemans et al., 2004). More, since the distance to some of the AM CVn stars 
is known with sufficient accuracy, they could be used as verification binaries for LISA (Stroeer and 
Vecchio, 2006). 

AM CVn stars, despite small number of identified and candidate objects (~ 20, see the list with 
parameters of the systems and references at 

www . astro . ru . nl/~nelemans/dokuwiki/doku . php?id=verif ication_binaries : am_cvn_stars) , 

are of great interest not only because they are potential sources of detectable gravitational waves, but 
also, for instance, because their formation involves evolution in common envelopes, their accretion 
1 For realization of this scenario angular momentum loss by magnetically coupled stellar wind is also necessary. 



disks may consist of helium or helium-carbon-oxygen mixture. AM CVn stars are potential pro- 
genitors of SN la (Tutukov and Yungelson, 1979; see also Solheim and Yungelson, 2005), of still 
hypothetical SN .la (Bildsten et al., 2007), and of other explosive phenomena associated with accu- 
mulation of helium at the surface of white dwarfs (see, e.g., Tutukov and Yungelson, 1981; Iben and 
Tutukov, 1991). Note also, that it is expected that ~ 100 most tight Galactic AM CVn stars might be 
observed both in gravitational waves and in the electromagnetic spectrum (Nelemans et al., 2004). 

It is still unclear which of the above-mentioned scenarios for formation of AM CVn stars acts in 
the Nature. All scenarios for formation of AM CVn stars imply that their progenitors pass through 
one or two common-envelope stages, but since clear understanding of the processes occurring in the 
latter is absent, the efficiency of common envelopes ejection remains a crucial parameter of scenarios 
that defines final separation of components and, hence, possibility of formation of a semidetached 
system. Also, conditions for stable mass transfer after RLOF by a white dwarf are not clear (see, 
e.g., Marsh et al., 2004; Gokhale et al., 2007; Motl et al., 2007). On the other hand, for instance, 
if AM CVn stars descend from strongly evolved hydrogen-rich cataclysmic variables, only a minor 
fraction of them evolves to P OI b < 25 min. which are typical for a considerable number of the 
AM CVn stars; more, hydrogen which has to be present in the spectra of accretion disks of most of 
such systems is not observed as yet. Nevertheless, it is possible that all formation channels contribute 
to the population of AM CVn stars. Resolution of the problem of the origin of AM CVn stars is also 
important because observational estimates of their Galactic population are by an order of magnitude 
below model estimates based on the above-described models (Roelofs et al., 2007a,b). This may point 
to certain flaws in the understanding of evolution of close binaries. Note, however, that the "deficit" of 
observed AM CVn stars may be due to numerous selection effects (see detailed discussion in Roelofs 
et al., 2007a,c and Anderson et al. 2005, 2008). 

Nelemans and Tout (2003) noticed that the chemistry of accretion disk in an AM CVn type star 
may serve as an identifier of its origin. If a helium white dwarf serves as the donor, the disk has to 
contain He and other H-burning products. In the case of helium- star donor, the disk, can, along to He, 
contain CNO-cycle and He-burning products. If the donor is an evolved main-sequence star, the disk 
may contain H and H-burning products. In the first case abundances depend on the mass of the white 
dwarf progenitor, while in other cases it depends also on how far was the star evolved prior to RLOF. 

Evolution of CB with low-mass He-donors was computed by Savonije et. al. (1986), Tutukov 
and Fedorova (1989), Ergma and Fedorova (1990). These papers were focused, mainly, on mass- 
transfer rates, ranges of orbital periods, since results of computations were applied to X-ray systems 
and SN la progenitors. In the present paper we carry out a systematic investigation of the evolution of 
helium stars in semidetached binaries as a function of their mass, evolutionary state at RLOF, and the 
parameters of a binary that contains helium star. Special attention is paid to the chemical composition 
of the matter lost by helium stars. In the forthcoming study (Nelemans and Yungelson, in prep.) these 
results will be applied for the analysis of the origin of observed AM CVn stars. 



2 Method of computations 



For our computations we have used a specially adapted version of P.P. Eggleton evolutionary 
code (1971; private comm. 2006). Equation of state, opacity tables, and other input data are de- 
scribed by Pols et al. (1995). Information on the latest modifications of the code may be found at 
www.ast.cam.ac.uk/stars/. In the code, abundances of ! H, 4 He, 12 C, 14 N, 16 0, 20 Ne, 24 Mg are 
computed. The nuclear reactions network is given in Table 1 of Pols et al. (1995). Nuclear reac- 
tions rates are taken after Caughlan and Fowler (1988), with exception of 12 C(a, y) 16 reaction for 
which the data of Caughlan et al. (1985) are used. Homogeneous helium models have mass frac- 
tions of helium Y=0.98, carbon X c =0.00019, nitrogen X N =0.01315, oxygen X Q = 0.00072, neon 
X Ne =0.00185, magnesium X Mg =0.00068. 

It was assumed that mass-exchange is conservative. Angular momentum loss via GWR was taken 
into account using standard Landau and Lifshitz (1971) formula 

'J\ 32 G 3 M He M wd (M He + M wd ) 
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Here M He and M wd are the masses of components, a is their separation. 
Mass loss rate M He is related to 7/7 as 
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where £(M He ) = dlnR/dlnMue- The term in brackets is ~ 1. 



3 Results 

3.1 Evolution of helium stars 

Following the model of the Galactic population of AM CVn stars (Nelemans et al., 2001), we 
have considered as typical progenitors of helium-donor AM CVn stars the systems with masses M He + 
M wd = (0.35 + 0.5) M Q , (0.4 + 0.6) M Q , and (0.4 + 0.8) M Q . In addition, we have considered a pair 
(0.65+0.8) M Q . If helium-accreting white dwarfs avoid double- (or "edge-lit-") detonation which may 
destroy the dwarf, such a system may belong to progenitors of AM CVn stars, but not to the most 
"fertile" of them. As the second parameter of computations we considered initial orbital period of the 
system, i.e., the period of the system immediately after completion of the common envelope stage that 
resulted in formation of a helium star. The range of P rb,o f° r every system was chosen in such a way 
that the set of computed models included both stars that filled their Roche lobe virtually unevolved 
and stars that had at the instant of RLOF radii close to the maximum of the radii of low-mass helium 
stars which are attained at Y c 0.1 (see Table. [Hand Fig.Q]). The maximum period in the set of initial 
orbital periods for every computed M He + M wd pair is, in fact, the limiting initial period which still 



Table 1. Parameters of calculated evolutionary tracks. Mj, M a - initial masses of donor and accretor, Po - 
initial orbital period, t c - time to RLOF , P c - orbital period at RLOF, Y c - central abundance of He at RLOF, 
tf - the age of the last computed model, Pf- final orbital period, Mjf - mass of the last computed model of the 
donor, Y s f -surface He-abundance of the last computed model. 
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0.35 


0.5 


80 


110.71 


18.14 


0.774 


492.51 


40.81 


0.027 


0.723 


5 


0.35 


0.5 


100 


202.32 


19.71 


0.642 


615.18 


40.35 
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0.6 


60 


50.78 


17.02 


0.871 


426.05 


41.18 


0.028 


0.825 


12 


0.40 


0.6 


80 


84.52 


22.42 
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0.495 


14 


0.40 
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0.021 
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0.057 


23 


0.65 


0.8 


35 


0.17 


34.54 


0.976 


324.85 


43.96 


0.027 


0.856 


24 


0.65 


0.8 


40 


2.06 


35.29 


0.928 


340.97 


44.04 


0.026 


0.806 


25 


0.65 
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allows RLOF and formation of an AM CVn star. Nuclear evolution of He-stars with mass M < 0.8 M ( 
terminates after helium burning stage and they evolve directly into white dwarfs (Paczyhski, 1971). 
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Figure 1. Dependence on time of abundances of He, C, N, O, and Ne in the center of 0.35 M , 0.40 M Q and 
0.65 M Q stars (top to bottom). The arrows mark the epochs of RLOF by stellar models computed in the paper 
(Table ffl). 

The lifetime of helium stars is short. According to our computations for (0.35 - 0.65) M Q mass 
range 

t Ue * 10 6 - 95 M£\l + Mf£ 4 ), (3) 

where time is in yr, masses in M Q . During He-burning, the radii of helium stars increase by ~ 30% 
only. Therefore, the crucial factor which defines the range of the orbital periods of the precursors of 
AM CVn stars is AML. In the Nelemans et. al. (2001) model, in 90% of the systems helium stars 
overflow their Roche lobes during first - 50% of their lifetime. Figure Q] shows that this time-span 
corresponds to the reduction of Y c to - 0.5. Thus, most of the AM CVn systems might have initial 
orbital periods up to 100 - 120 min. 

Evolution of mass-losing stars is followed to M (0.02 - 0.03) M Q . Continuation of computa- 
tions is hampered by the absence of adequate EOS and low-temperature opacity tables in our code. 
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Figure 2. Dependence of mass-loss rate on the period for computed evolutionary sequences for the systems 
with initial mass of components 0.35 M Q + 0.5 M Q , 0.40 M Q + 0.8 M Q , 0.65 M G + 0.8 M Q (top to bottom). 

However, for the analysis of the scenarios of formation and evolution of AM CVn stars which involves 
information on the elemental abundances in the donor-star this is not important, since the chemical 
composition of the stars in our case becomes "frozen" even before the minimum of P orb is reached, 
because of the drop of the temperature in the core and switch-off of nuclear burning (see below). 

Figure |2] shows dependence of the mass-loss rate by Roche-lobes filling helium stars on P orb . 
The behavior of low-mass components of CB evolving under the influence of AML via GWR, like 
for "usual" hydrogen-rich cataclysmic variables, is defined by the relation between the timescale of 
stellar evolution r ev , thermal timescale of the star tkh and AML timescale t G w (Faulkner, 1971; 
Paczyhski, 1981; Savonije et al., 1986). The relation between the timescales itself depends on the 
mass of the star, its evolutionary state at RLOF and total mass of the binary (see Fig.©. 

In the system (0.35+0.5) M , P orb0 =20min. with initially virtually unevolved donor, tkh > 1gw 
after the RLOF. Due to AML separation of components and Roche-lobe radius decrease. The star 
which is not in thermal equilibrium reacts to the mass loss by an increase of mass-loss rate. But since 
mass-transfer acts in opposite direction (moves components apart), M does not change significantly. 
A considerable fraction of the energy generated by nuclear burning is absorbed in the envelope of 
the star, resulting in decrease of luminosity and further growth of tkh- With decrease of the stellar 
mass the significance of nuclear burning is rapidly diminishing, convective core disappears, while 




P, min P, min 

Figure 3. Variation of thermal timescale of helium star tkh, angular momentum loss timescale tgw and donor 
radius (in R Q ) in the binaries with initial parameters M2,o - 0.35 M Q , M\$ - 0.5 M Q and P O rb,0 = 20 min. (left 
panel) and M2,o = 0.65 M Q , M\fl - 0.8 M G and P m b,o = 60 min. (right panel). A constant is added to the radii 
for combination of the curves in the same plot. 

a surface convective zone appears; during this stage t K h » ^gw- As P orb drops, r G w continues to 
decrease, resulting in increase of M. At certain instant the effect of mass-transfer starts to dominate, 
orbital period reaches the minimum, while M - the maximum. With increase of P orb the timescale of 
AML increases, but since the star becomes more degenerate and has convective envelope, the power 
of M - R relation is negative and mass loss continues, but M rapidly drops. 

If the donor is initially more massive and more evolved at the instant of RLOF, it retains thermal 
equilibrium for a longer time and M is initially completely defined by AML. But trh increases with 
decrease of M He , while tgw decreases and, starting from a certain moment, evolution proceeds like in 
the lower-mass and less-evolved system described above. Evolution of characteristic timescales for 
the system (0.65+0.8) M , i> or b,o =60 min. is similar to that described by Savonije et al. (1986) for a 
(0.6+1.4) M system with Y c = 0.28 at RLOF which is usually quoted as "typical". 

The minimum periods of the systems with helium donors are confined to a narrow range of 
^"orb = (9.3 — 10. 9) min, masses of stars at /Vb.mm are 0.20 M Q - 0.26 M Q . The models that overfilled 
their Roche lobes in the advanced stages of evolution reach lower periods and have at minimum period 
lower masses than initially less evolved donors. This may be related to the lower abundance of He 
and enhanced abundance of heavy elements. 

3.2 Mass-radius relation 

All existing computations of the evolution of initially non-degenerate helium donors are discon- 
tinued at M He - 0.02 M Q . However, starting from P orb ^ 20 min. these stars are, in fact, weakly 
degenerate homogeneous objects (the parameter of degeneracy \p « 8.33 x 10 7 p c T c T 3/2 ~ 10). This 
allows to consider qualitatively their further evolution using results of computations for mass-losing 
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Figure 4. Dependence of model radii on stellar mass for the systems 0.35 M Q +0.5 M Q (upper panel), 0.40 M Q + 
0.6 M Q (middle panel), 0.65 M Q + 0.8 M Q (lower panel). Less evolved systems are at the top of every panel. In 
the upper panel upper dashed line shows M-R relation dU, lower dashed line shows M-R relation for 0.3 M Q , 
log \]/q - 1.1 helium white dwarf (Deloye et al., 2007) and dash-dotted line shows the fit to M - R relation for 
zero-temperature He white dwarfs suggested in the latter paper. 

arbitrary degenerate helium dwarfs (Deloye et al., 2007). Deloye and coauthors have shown that, 
when the mass of the dwarf decreases to - 0.01 - 0.03 M , their tkh becomes comparable to the 
timescale of mass loss r m and further tkh < T m - The stage of adiabatic expansion terminates, the 
donor gradually cools, becomes more degenerate and M-R relation approaches relation for fully 
degenerate configurations. Expected picture of evolution is confirmed by Fig. 0] which shows vari- 
ation of stellar radii with decrease of mass. In addition to our data, we show in Fig. |4] mass-radius 
relation for degenerate helium white dwarf with initial mass of 0.3 M Q and initial degeneracy param- 
eter logi/r = 1.1 (the least degenerate initially model, C. Deloye priv. comm.) and M-R relation 
for zero-temperature He white dwarfs (Deloye et al. 2007). It is seen that the morphology of M - R 
curves for helium stars and helium dwarfs is similar and, hence, one may expect that the remnants 
of helium stars will gradually approach the curve for the remnants of helium dwarfs^. Note, Deloye 
et al. in fact predicted this effect, based on consideration of the physics of cooling objects and they 

2 When comparing our results with those of Deloye et. al. one must have in mind the differences in initial masses of 
models and degeneracy parameters, as well as in the input parameters of different codes. 
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Figure 5. Dependence of P rt>/A)rb on P or b for the 0.35 M G + 0.5 M G , P O rt>,0 - 20 min. system. Solid line 
- model calculations, upper dashed line - dependence obtained using M - R relation ©, lower dashed line - 
dependence obtained by means of relation ©. 

noticed also that some of the evolutionary tracks published by Tutukov and Fedorova (1989) have 
features described above. 

Our results, combined with results of Deloye et al. (2007) suggest that at P orb ~ 40 min. two 
families of AM CVn stars - the "white dwarf" one and the "helium star" one - merge and the origin 
of the donors may be then identified by their chemical composition only (see §@]). 

Our conclusions concerning M - R relation may have certain implications for population models 
of AM CVn stars. Existing models (Tutukov and Yungelson, 1996; Hils and Bender 2000; Nelemans 
et al., 2001, 2004) were carried out under assumption that the power of M - R relation is constant, 
keeps also for M He ~ 0.02 M Q and that the evolution of stars is restricted by the Hubble time only. 
This assumption allowed to use for modeling Eq. © with constant £(M He ). Such an approach is, ap- 
parently, not justified and one has to use an M- R relation that takes into account growing degeneracy 
of donors with mass loss. 

Figure |4] shows that the M - R relation (in solar units) 

R « 10" 1367 M H ^ 062 , (4) 

suggested by Nelemans et al. (2001) using results of Tutukov and Fedorova (1989) does not agree 
well with the results of more modern computations and is not valid for P orb > 35 min. 

In the period range P mb w 10-35 min which hosts 10 out of 17 AM CVn stars with estimated 
orbital periods, the radii of the models with initial masses (0.35-0.40) M Q may be approximated as 

/ p \-0-05 , n q <- v 0.345 
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Figure 6. Orbital period - donor mass relations. For every given combination of initial masses of components 
we show P or b - M relation for the system with initially least evolved and most evolved donor. 0.35 M Q +0.5 M 
- thin solid lines, 0.40 M Q +0.6 M Q - dashed lines, 0.65 M Q +0.8 M Q - dotted lines. For initially evolved donors 
the lines practically coincide. Dash-dot-dot-dot line shows P 0Tb -M relation for 0.35 M Q +0.5 M Q , P O rt>,0 system 
based on approximation ©. Dot-dash line shows P or b - M relation for Mq = 0.3 M Q , log ifro - 1 . 1 helium white 
dwarf from Deloye et al. (2007). Vertical bars show the ranges of donor-mass estimates for AM CVn, HP Lib, 
CR Boo, V803 Cen, SDSS J0926+3624, and GP Com (see references in the text). 

where P rbfl and M He ,o are initial orbital period of the system and initial mass of the donor and M and 
R are in solar units, P 0I \, - in min. 

The time spent by a star in the certain range of orbital periods is proportional to /Vb/^Vb- In 
Fig. |5] we compare dependencies of P or b/P rb on P or \, for the system (0.35+0.50) M Q , P rb,o = 20 
obtained in evolutionary computations and obtained by means of Eqs. © and ©. It is evident that 
in the P mb « (10 - 35) min. range Eq. © results in the 20-25% difference in the number of systems. 
Note, however, that this discrepancy is not very significant if one takes into account all uncertainties 
involved in the modeling of the population of AM CVn stars. 

3.3 Masses of donors in AM CVn stars 

Relatively accurate masses of components are derived only for SDSS J0926+3624 with P orb =28.3 min, 
a unique partially eclipsing AM CVn type system (Marsh et al., 2007): M wd = 0.84 + 0.05 M Q , 
M He = 0.029 + 0.002 Mq. In the P - M diagram (Fig. [6]) the donor of SDSS J0926+3624 is located 
both below the curve corresponding to the initially most evolved nondegenerate helium donors and the 
curve that describes the least degenerate white dwarfs. Apparently, this system belongs to the latter 
family of stars (see also Deloye et al., 2007). Note however, that at the moment when P orb =2S.3 min, 




the remnants of initially substantially evolved donors become oxygen-neon white dwarfs (see Figs. [8] 
— ITOl). Though existing models of the population of AM CVn stars suggest that the probability of 
RLOF by a far-evolved helium star is low, only analysis of chemical composition of the donor will be 
able to identify the scenario of the origin of SDSS J0926+3624. 

We show in Fig. [6] also mass estimates for several other AM CVn stars after Roelofs et al. 
(2006; 2007b). For AM CVn the masses of components are based on the kinematical data: M W( j = 
(0.68 + 0.06) M Q , M He = 0.125 + 0.012 M . For HP Lib, CR Boo, and V803 Cen masses of the donors 
are derived indirectly, assuming that for these systems one may use the relation between superhump 
period excess and mass ratio of components derived by Roelofs et al. (2006) for AM CVn: s(q) = 
0.12^. For GP Com the range of possible donor mass is based on the limits for its luminosity and 
assumptions about the nature of the star. The estimated masses of the donors agree with assumption 
that they descend from nondegenerate helium stars (see also Roelofs et al., 2007a). Note however, 
that there are no traces of He-burning products in the spectra of these stars. 

Spectral lines of nitrogen, neon, and helium were discovered in GP Com (Lambert and Slovak, 
1981; Marsh et al., 1991, 1995; Strohmayer, 2004). In particular, Strohmayer found the following 
abundances: Y = 0.977 + 0.002, X N = 1.7 + 0.1 x 10~ 2 , X Q = 2.2 + 0.3 x 10 3 , X Ne = 3.7 + 0.2 x 10 3 , 
X c < 2 x 10~ 3 , X Mg < 1.6 x 10~ 4 . Considerable excess of nitrogen compared to oxygen and carbon 
points to the enrichment of stellar matter by the products of CNO-cycle. Such abundances may be 
typical for the remnants of the least massive possible progenitors (M He ,o ~ 0.3 M ) which overfilled 
Roche lobe soon after TAMS, before helium started to burn in their interiors. In such stars abundances 
in the matter lost by them is virtually the same during whole course of evolution and close to the ones 
found by Strohmayer. 

3.4 The system (0.65+0.8) M , P orb () = 90 min. 

The system with initial parameters (0.65+0.8) M Q , P O rb,0 = 90 min. is of special interest because 
of its "nonstandard" evolution which may illustrate transitional behavior between "white dwarf" and 
"helium star" scenarios of formation of AM CVn stars. At the instant of RLOF abundances in the 
center of the donor are Y « 0.19, X c ~ 0.45, X « 0.34, nitrogen is destroyed already. At difference 
to lower mass stars in which nuclear burning terminates after loss of several hundredth of M Q , in 
this more massive and hotter star it lasts longer. When the mass decreases to M 2 « 0.52M o and Y c 
becomes « 0.008, the star begins to contract. Mass-exchange interrupts for 2.2 Myr and resumes in 
the helium-shell burning stage (Fig. [7]). After the loss of additional ^0.1 M Q nuclear burning ceases 
and the star is now a "hybrid" white dwarf with =* 0.3 M Q carbon-oxygen-neon core and ^0.1 M Q 
helium envelope with an admixture of C and O. The timescale of AML r G w continues to decrease, 
while t K h increases. Evolution of the donor is accompanied by its cooling and decrease of radius. 
The dwarf experiences a stage of adiabatic contraction similar to the one that is typical for the initial 
stages of mass-exchange in the systems with degenerate donors (Deloye et al., 2007). In the last 
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Figure 7. Lower panel: dependence of mass, radius and central helium abundance on time for the model with 
initial mass Mif) — 0.65 Mq in the system with M\fi — 0.8 Mq, ^orKO — 90 mm. Upper panel; Variations of 
angular momentum loss timescale and thermal timescale of the donor. 

computed model of the sequence, P orb =4.3 min, M He =0.375 M Q , degeneracy parameter of the center 
of the star is i/r « 14. One may expect that further evolution will be qualitatively similar to the 
evolution of helium white dwarf-donors: contraction phase will be followed by expansion which will 
continue until mass will decrease to 0.01 M Q and thermal timescale will become comparable to the 
mas-loss timescale and a stage of contraction to the fully degenerate configuration will ensue. The 
minimum P orb of the system will be several minutes. Note, that at difference to He white-dwarf donors 
which retain their original chemical composition during evolution, in the case of He-star descendant 
under consideration, the donor will become an oxygen-carbon white dwarf, with an admixture of neon 
(X « 0.27, X c ~ 0.71, see Fig.Q). 

Total mass of this system is 1.45M . If unstable He-burning at the surface of white dwarf does 
not decrease its mass below Chandrasekhar limiting mass, such a system may be a SN la progenitor. 

4 Evolution of stellar chemical composition 

Figures [8] -[Kj] show the evolution of chemical composition of matter lost by stars upon RLOF 
As we mentioned before, in the stars with initial masses 0.35 and 0.40 M nuclear burning termi- 
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Figure 8. Variations of He, C, N, and abundances in the accreted matter in the system with initial masses of 
components 0.35 M Q and 0.5 M Q . Solid lines are for systems with initial orbital periods P ib,o =20, 40 and 60 
min., dashed lines - for -P O rb,0 =80 min., dash-dot lines - for P O rb,0 =100 min., dotted line - for P O rb,0 =120min., 
dash-dot-dot-dot line - for P O rb,0 =140 and 144 min. For the system with P O rb,0 =20 min. He-abundance 
virtually does not change. All abundances do not change any more at /> or b ^25 min. 

nates soon after RLOF. Simultaneously, disappears the convective core. When P orb is still larger 
than P tb,min> abundances in the lost matter are virtually constant and correspond to the abundances 
in the matter that experienced hydrogen burning via CNO-cycle. In the period range from P rb,mm to 
Port. » (15 - 20) min. the convective core of initial model is uncovered. The donors in initially most 
wide systems are an exception, since they had large cores and chemical composition of the lost matter 
starts to change already before P rb, m m- But this change happens over such a short time that it is highly 
unlikely to observe a star in this "transition" state. Thus, in the range of P orb that harbors almost all 
known AM CVn stars chemical composition of the matter accreted by white dwarf is defined by the 
extent of helium exhaustion at the instant of RLOF (compare Fig. Q]and Figs. [81— ITDl). 

As we mentioned above, it is expected that the overwhelming majority of precursors of AM CVn 
stars in the "helium star" channel of formation were formed with P mb < 100 - 120 min. Then, 
typical abundances in the transferred matter are Y > 0.4, 2 x 10 -4 < X c ~ 0.3, 7 x 10" 4 < X < 0.25, 
5 x 10" 4 ~ X N < 0.5 x 10" 2 . However, one cannot exclude that abundance of nitrogen is much lower 
or it is absent at all; diminished (or zero) X N has to correlate with relatively low Y and enhanced X c 
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Figure 9. Same as in Fig.[8]for the system with initial masses of components 0.4 M Q and 0.6 M Q . Solid lines 
show abundances for the systems with P O rt>,0 =20, 40 and 60 min., dashed lines - for P rb,o = 80 min., dash-dot 
lines - for P O rb,0 =100 min., dashed line - for P O rb,0 =120 min., dash-dot-dot-dot line - for P m b,o =130 min. 

and X Q . 

We did not plot in Figs. [8]-[l0] variations of the Ne abundance. It also changes rapidly by P OT b, min 
and at P orb >20 min. neon abundance 2 x 1(T 3 < X Ne < 2 x 10" 2 . Variation of X Ne is, evidently, too 
small to serve as identifier of the scenarios of AM CVn stars formation. 

If the donor in an AM CVn type system was initially a white dwarf, there are no reasons to expect 
that abundances in the matter lost by it vary in the course of evolution. They must to correspond to 
the "standard" chemical composition of a stellar core that experienced hydrogen burning, with small 
variations due to the differences in the masses of precursors of white dwarfs. Similar abundances will 
have matter lost by He-stars that filled Roche lobes soon after TAMS. However, since initial Y « 1, 
while initial X c ~ 10" 4 , Z N ~ 10~ 2 , Xo ~ 10" 3 , detection of even slight enrichment in C and O or 
impoverishment in N may indicate possibility of formation of AM CVn stars with He-star donors. 

5 Conclusion 

Above, we presented results of evolutionary computations for semidetached low-mass helium 
components in the systems with white dwarf accretors. Evolution was considered as conservative 
in mass, but accompanied by angular momentum loss via gravitational waves radiation. We may 
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Figure 10. Same as in Fig.[8]for the system with initial masses of components 0.65 M Q and 0.8 M Q . Solid lines 
show abundances for the systems with P O rb,0 -35 and 40 min., dashed lines - for P rb,o =60 min., dash-dot 
lines - for P O rb,0 =80 min., dashed line - for P O rb,0 =85 min., dash-dot-dot-dot line - for P O rb,0 =90 min. 



summarize our main results as follows. 

Evolution of binaries under consideration only weakly depends on the mass of helium star, total 
mass of the system, and evolutionary state of the donor at the instant of RLOF. The minimum P 0I \,= 
9.3 - 10.9 min, the masses of the donor at P rb,min are 0.20 M ~ 0.26 M Q (with exception of initially 
very far-evolved system (0.65 + 0.8) M Q , P or b,o = 90 min.) which, in fact, represents a transitional case 
between systems with initially nondegenerate and degenerate donors). In the period range 20 - 40 
min, which may be adequately described by our models and which hosts a significant fraction of all 
observed AM CVn stars, the scatter of M for a given P orb does not exceed a factor ~ 2.5. 

Our computations are limited by M He > 0.02 M Q . With further decrease of stellar mass thermal 
timescale of the donor becomes shorter than the angular momentum loss timescale, the stage of adi- 
abatic donor expansion comes to the end, and the matter of the donor becomes more degenerate as it 
cools down. Qualitative comparison with the computations for initially arbitrary degenerate helium 
white dwarfs allows to suggest that at P orb > (40-45) min "white-dwarf" and "helium-star" families 
of AM CVn stars merge and only analysis of the chemical composition of the matter lost by the donor 
is able to discriminate between different scenarios of formation. Such an analysis is carried out in 
continuation of the present study. This result also shows that existing theoretical models of the pop- 



ulation of AM CVn stars might need certain revision, since they were carried out under assumption 
that helium stars evolve till Hubble time with the same M - R relation. 

Initial periods of the most typical, according to current understanding, progenitors of AM CVn 
stars with initially nondegenerate helium donors are confined to about (20 - 120) min. This implies 
a very narrow range of initial separations of components (for instance, 2 R Q - 7 R Q for M He + M wd = 
0.35 M Q + 0.5 M Q pair). This interval is defined by a combination of different parameters. A binary 
passes through two common envelope stages. The "standard" equation for the variation of orbital 
separation of components based on energy balance between binding energy of the mass-losing star 
and orbital energy of the system (Webbink, 1984; de Kool et al., 1987) is 

a f Mi jC 
at ~ Mi 

where the indexes i and / label initial and final separation of components, a - is the parameter of 
common envelope efficiency, A - is the parameter of the binding energy of the stellar envelope, M\ 
and M l c are initial mass of mass-losing star and the mass of its remnant, r l L is dimensionless radius 
of the star at the beginning of mass transfer, M 2 is the mass of companion. 

For every common-envelope stage the values of a and A are, most probably, different and depend 
on a, at the beginning of the stage; the value of cif after the second common envelope depends on 
"initial-final" mass relations for progenitors of white dwarf and helium star, which, in their own turn, 
depend on their evolutionary state at the instant of RLOF. This suggests that formation of AM CVn 
stars needs a very fine "tuning" of evolutionary parameters. The treatment of stellar evolution in- 
corporated in the population synthesis codes that predict existence of AM CVn stars with initially 
nondegenerate He-donors is, at the moment, inevitably too crude for account of all subtleties of pa- 
rameters. If AM CVn stars with confidently established presence of He-burning products signatures 
in their spectra will not be detected, this may mean that the "necessary" combination of parameters 
is not realized in the Nature or this combination is such, that all He-donors overflow Roche lobes 
prior to ignition of He or when He burns very weakly. Such a possibility is suggested by positions 
of AM CVn and V803 Cen in the period-mass diagram (Fig. |6]). Since initial abundance of He in 
stars is close to 1, while abundances of C, O, N are by 2 to 4 orders of magnitude lower, even slight 
enhancements of X c and X or reduction of X N may indicate possibility of formation of AM CVn 
systems with He-star donors. 

Above-discussed scenario of evolution needs some caveats. Evolution of every considered sys- 
tem has two phases - prior and after period minimum. In the systems with M 2j o ~ (0.35 - 0.40) 
and Porb.o - 120 min. which we consider as typical, in the first stage of mass-exchange 2.5 x 
10~ 8 < M < 10~ 7 M Q yr 1 . Livne (1990), Livne and Glasner (1991), Livne and Arnett (1995) and 
other authors have shown that for accretion rate M ~ 10~ 8 M G yr' 1 accumulation of ~ 0.1 M© of 
He at the surface of CO white dwarf may result in detonation of He which may initiate detonation 
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of carbon in the center of white dwarf and, presumably, SN Ia_j. Modeling by means of population 
synthesis has shown that ELDs may be dominating mechanism of SN la in young (< 1 Gyr) popu- 
lations (see, e.g., Branch et al., 1995; Hurley et al., 2002; Yungelson, 2005b). However, models of 
light-curves (Hoeflich and Khokhlov, 1996; Hoeflich et al., 1996) and spectra (Nugent et al., 1997) 
of objects experiencing ELDs revealed that they disagree with observations. But this conclusion is 
by no means final in the absence of more detailed calculations. Double-detonation, if it occurs, may 
prevent transformation into AM CVn stars of some (~ 40%) of their potential precursors (Nelemans 
et al., 2001). Note, for instance, that double-detonation may prevent formation of an AM CVn star 
with initial parameters M He + M wd = (0.65+0.80) M Q . It is possible that instead of a detonation of a 
massive He-layer recurrent Nova-scale helium flashes occur (see, e.g., Bildsten et al., 2007 and ref- 
erences therein). The flashes may be accompanied by mass and momentum loss from the system and 
change the behavior of evolutionary sequences. The problem of unstable He-burning at the surface 
of accreting white dwarfs and its possible impact on the evolution of AM CVn stars needs additional 
investigation. 

Availability of results: Detailed results of computations may be found at 
www . inasan . ru/~lry/HELIUM_STARS/. 
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